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ABSTRACT 

We evaluate the effects of high-energy ionizing particles on the Si:As impurity band 
conduction (IBC) mid-infrared detector on board AKARI, the Japanese infrared as- 
tronomical satellite. IBC-type detectors are known to be little influenced by ionizing 
radiation. However we find that the detector is significantly affected by in-orbit ionizing 
radiation even after spikes induced by ionizing particles are removed. The effects are 
described as changes mostly in the offset of detector output, but not in the gain. We 
conclude that the changes in the offset are caused mainly by increase in dark current. 
We establish a method to correct these ionizing radiation effects. The method is essen- 
tial to improve the quality and to increase the sky coverage of the A KARI mid-infrared 
all-sky-survey map. 



1. Introduction 



The effects of high-energy ionizing particles (hereafter we simply call them ionizing radia- 
tion effects) on infrared detectors of extrinsic semiconductors are serious for infrared astronomical 
observations in space. The ionizing radiation effects cause changes in detector sensitivity and after- 
effects which can last for a certain amount of time after the particle h it. For infrared astronomical 



satellites such a s Infrared Space Observa tory (ISO; iKessler et al J 1 19961 ) and Infrared Astoronomical 



Moselev et al. 


1989; 


Nieminen 


2003; 


Her as et al. 


2003) 



effects from their observational data (e.g. 
This problem is particularly serious in the south Atlantic anomaly (SAA). The SAA has the local 
depression in the Earth magnetic field intensity over Brazil, where many high-energy protons and 
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electrons are trapped even at low sky altitudes. When a satellite passes through the SAA, a great 
amount of ionizing radiation is irradiated onto detectors. 

AKARI, the first Japanese infrared astronomical satellite, was launched on February 21, 2006 
(UT) into a su n-synchronous polar or bit at an altitude of 700 km with a period of about 100 minutes 
per revolution ([Murakami et al.ll2007l ). AKARI made an all-sky survey in the six photometric bands 
centered at the wavelengths of 9 and 18 /jm with the Infrared Camera (IRC; Wada et al. 2003; 
Onaka et al. 2007) and 65, 90, 140, and 160 fim with the Far-Infrared Surveyor (FIS; Kawada 



et al . 2007) within a peri od of almost one and a half years (jKataza et al.l |2010| ; lYamamura et al 



2010i : llshihara et al.ll2010i ). The all-sky survey was performed during the lifetime of liquid helium 
cryogen between May 8, 2006 and August 28, 2007. As a result, A KARI surveyed more than 96 
and 97 % of the whole sky twice or more in the S9W and L18W bands, respectively. During the 
survey, AKARI passed through the SAA three or four times a day. 

For the far-infrared survey, the AKARI/FIS adopted Ge:Ga photoconductive devices, which 
made the all-sky surveys in the 65, 90, 140, and 160 /im bands. The Ge:Ga detectors were seriously 
influenced by the ionizing radiatio n, which caused undesirable changes in detector responsivity and 
degradation in noise performance (jKaneda et al.ll2006l : ISuzuki et al.ll2008l ). Thus a curing operation 
was performed every time the satellite passed through the SAA. The FIS did not carry out pointed 
observations during the passage of the SAA. 

The mid-infrared channels of the AKARI/IRC, which made the all-sky surveys at th e 9 //m and 



18 /xm bands, adopted the Si:As/CRC-744 infrared arrays manufactured by Raytheon (|Lum et al 



1993 



Estrada et al. 


1998; 


Ando et al. 


2003). 


(Szmulowicz & Madarsz 


1987 


) , which is 



like the above non-IBC-type Ge:Ga detectors ([Gordon et al.ll2005l ). However, spikes are inevitably 
induced by ionizing particle hitting events in the detector output, which resemble the detection of 
astronomical poin t sources. Therefore , algor ithm was developed to remove the spikes induced by 



ionizing particles (jlshihara et al 



2006 



20101 ). So far, it has not been quantitatively investigated 
whether or not there are any significant changes in detector performance other than the spikes due 
to in-orbit ionizing radiation. 

In this paper, we evaluate ionizing radiation effects on the Si:As IBC mid-infrared detector 
on board AKARI. We find that the detector is significantly affected by in-orbit ionizing radiation 
even after spikes induced by ionizing particles are removed. The effects are described as changes 
mostly in the offset of detector output, but not in the gain, where the offset value increases with 
the particle hitting rate. We establish a method of correction, by which we can utilize the data 
affected by the ionizing radiation to make diffuse sky maps. 
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Fig. 1— AKARI 18 /an maps of the 1°.7 x 1°.7 area centered at (R.A., Dec.)=(357°.7, -19°. 7): 
(a) original image created from data including those affected by the ionizing radiation inside the 
SAA, (b) the image created by not using the SAA-affected data, and (c) the image created after 
correcting the SAA-affected data by using the correction method presented in this paper. The red 
boxs in the panel (a) show the image units (see text) with particle hitting rates higher than 5 events 
unit -1 , which are defined as SAA-affected data. 



2. Data reduction 



We use the all-sky survey data in the 9 /im and 18 /im bands to investigate ionizing radiation 
effects on the Si:As IBC detectors. The contents of this paper are based on the 18 /jm data; for the 
9 /im data, similar results are obtained and the same technique as described below can be applied. 

In the all-sky survey observations, two out of 256 rows in the 256 x 256 dete ctor array format 
are operated in a continuous non-destructive readout mode (jlshihara et al.ll2006l ). The pixel scales 
are 2'.'34 x 2'.'34 for the 9 /im band and 2'.'51 x 2 / .'39 for the 18 /im band. The sampling rate is set 
to be 22.72 Hz (one sampling per 4.4xl0 -2 s), while the scan speed of the satellite in the survey 
observation mode is 216" s" 1 . The effective pixel size for a single scan by a single row is 9 / .'36 x 9'.'36 
for the 9 /im band and 10'.' 4 x 9 / .'36 for the 18 /im band. The in-scan pixel size is fixed by the 
sampling rate, while the outputs of four neighboring pixels are coadded to meet the down-link rate 
requirements and still make effective observations. Finer spatial resolutions comparable to the point 
spread functions (P SFs) of 5'.'5 and 5'/ 7 in FWHM for the 9 /im and the 18 /im band, respectively 
(lOnaka et al.1 120071). are o btained in the data processing by combining images produced by the two 
rows (jlshihara et al.ll2010l ). 



All the pixels were reset at a rate of 0.074 Hz (one reset per 306 samplings) to discharge 
the photo-current integrated in detector capacitance. The region of 10' (= 2 / /34 x 256 pixels) 
x48'(=216"/0.074 Hz) was surveyed in this reset interval. Here and hereafter we define this region 
as an image unit. After we process the data with the pipeline used in generating the point source 
catalog, we calculate averaged pixel intensities per image unit. The spikes induced by ionizing 
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particle hitting are detected by millisecond con firmation by using the redundancy of the 2 rows, 



and removed from the data (jlshihara et al.ll2010l ). We define the number of the spikes detected per 



image unit as particle hitting rate per image unit. 

Figure [1] presents an example of AKARI 18 fim maps of the 1° .7 x 1° .7 area centered at (R.A., 
Dec.)=(357°.7, —19°. 7), where we have already removed the spikes induced by ionizing particles 
from the data. The size of the image bin is set to be 6'.'24 x 6'.'24. Figure is the image created 
by using the data including those taken during the passage of the SAA, while Fig. []J> is the image 
created by not using such data. The bright stripes are seen in the original image (Fig.[TJi). However, 
after the data taken during the passage of the SAA are excluded, the stripes disappear from the 
image (Fig. QJ>). We thus find that the ionizing radiation has significant effects on the IBC detector 
output even after the spikes are removed. Here and hereafter, the SAA-affected data are defined 
as those included in image units with particle hitting rates higher than 5 events unit" 1 . 

As shown in Fig.Qja, the ionizing radiation apparently increases the output level of the detector, 
which can be explained by change in either offset or gain of the output. Because A KARI observed 
the same region twice or more and the IRC was continuously operated around the SAA, we can 
compare SAA-affected and non-affected intensities to quantitatively evaluate the change by the 
ionizing radiation by using two datasets taken at the same sky position in different seasons. 

However, we cannot simply compare the data in different seasons because the Zodiacal light 
emission, a dominant foreground emission component in the mid-IR wavelength region, is known 
to show seasonal variations due to asymmetr y in distribution of Zodiacal dust with respect to the 



ecliptic plane and the circum-solar dust ring (jKelsall et al.lll998l ). Thus, we first develop a method 



to evaluate the seasonal variation of the foreground emission. 

The difference of the foreground emission (A/ z ) between two seasons due to the variation of 
the Zodiacal light is represented as, 



AI Z = h- I 2 , (1) 

where I\ and I2 are the averaged pixel intensities of image units at the same sky position observed 
in the seasons 1 and 2, respectively. Both I\ and I2 are not affected by the SAA. Figure [2] shows 
an example of AI Z as a function of the orbital phase of AKARI (the orbital period of AKARI is 
approximately 100 minutes). The time corresponds to the orbital phase, with which the terrestrial 
latitude of the satellite position and the ecliptic latitude of the observed sky change. The I\ values 
are taken from data obtained on June 10, 2006, while the I2 values are taken from data obtained 
half a year later, on December 10, 2006. Because the orbit of AKARI is sun-synchronous, almost 
the same sky positions are covered in these two observations. The plots A-E in Fig. [2] correspond 
to data in 5 sequential orbits. The maximum distance of sky positions at the same orbital phase 
between two adjacent orbits is 4 arcmin, and that between orbits within 24 hours is one degree. 
The plots B, C, and D correspond to data in the orbits passing the SAA. Large discrepancies 
between the estimated and the observed AI Z are seen during the passage of the SAA, where the 
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data cannot be used in estimating AI Z . They are due to the ionizing radiation effects. By assuming 
that AI Z has a smooth spatial distribution within one degree and a small time variation within 24 
hours, non-SAA-affected levels of AI Z are estimated by interpolating from the neighboring orbit 
data along the cross-scan direction. We do not use data when the resulting AI Z shows a time 
variation faster than 0.44 ADU per minute. By using these plots, we can distinguish the effects of 
ionizing radiation from the seasonal variation of the foreground emission. 



3. Results 

3.1. Behavior of ionizing radiation effects 

As shown in the previous section, the output levels of the detector increase systematically in 
the SAA even after the spikes are removed. Therefore we examine how they increase, i.e. by change 
in offset or gain of the output. 

The fractional increase of the output level due to the SAA, f{lf AA )^ is derived from the 
equation, 

ti jSAA-, _ {Ii AA ~ h ~ AJ Z ) 

J(h ) ~ JsAA ' 

where lf AA is the averaged pixel intensity of an image unit observed in the season 1 during the 
passage of the SAA. If the ionizing radiation effect is described by the change in the gain (1+a), 
f{lf AA ) is reduced to, 

= (1+ °»Vf- AJ - = (3) 
(1 + a)Ii 1 + a 

The scatter plot of f{lf AA ) versus lf AA should follow the line parallel to a horizontal axis, if a 
is constant. On the other hand, if the ionizing radiation effect is described by the offset the 
scatter plot follows a hyperbolic function as, 

SAA _ (h + (3)-I 2 - AI Z P 

Therefore, we can judge whether the ionizing radiation effect is the gain or the offset from the 
behavior of f{lf AA ) with respect to lf AA . 

Figure E] shows the plot of f{lf AA ) versus /^ AA for the three different ranges of particle hitting 
rates. Because these plots are better fitted by hyperbolic curves than by lines parallel to the 
horizontal axis, we conclude that the ionizing radiation effect is explained mostly by the change in 
the offset of detector output, but not in the gain. 

Figure U] compares the distribution map of the offset value of the detector output ((3=I SAA — I) 
with the contours of the particle hitting rate in the terrestrial coordinates. The distribution map of 
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the offset is spatially correlated well with the contour map of the particle hitting rate, and does not 
elongate along the satellite scan direction, indicating that after-effects due to the ionizing radiation 
are negligible, at least, on time scales longer than a few minutes. 



3.2. Correction of ionizing radiation effects 

We evaluate the ionizing radiation effects quantitatively to formulate a correction method. 
t igure plots offset values (/3=I SAA -I) as a function of the particle hitting rate. The plot shows a 
linear relation between the offset values and particle hitting rates. We obtain the linear coefficient 
of 0.191 ± 0.001 by a line fitting. 

Thus in order to make diffuse maps, we first obtain a particle hitting rate every image unit 
to be used. Then as a function of the particle hitting rate, we derive the corresponding offset 
values from the line in Fig. [5j and subtracted them from the image units affected by the SAA. As 
a result, we derive such a clean image as shown in Fig. [lb, which demonstrates the effectiveness 
of this correction method. We also apply this method to another region near the Galactic center 
(Fig. [6]). As can be confirmed from the figure, this method works out well in making diffuse maps 
with brighter intensities and larger dynamic ranges of brightness. Thus by making the most use of 
the data affected by the ionizing radiation with this method, we can increase the sky coverage and 
the quality of the AKARI mid-IR all-sky-survey map, especially the Galactic plane map. 



4. Discussion 



What causes the offset of the output level during the passage of the SAA? Possible scenarios are: 
(1) an increase in the number of very small spikes below the detection level, (2) an increase in dark 
current, or (3) a drift of the offset level of readout electronics. The last one is due to the change of the 
gate-sour ce voltage of a source- follower in the readout integrated circuit (ROIC) hy bridized to the 



detec tor (jWu et al.lll997l ). which is easily caused by the drift of an array temperature (jlshihara et al 



2003). To distinguish them, we compare the distribution of the pixel intensities of the image unit 
taken inside the SAA with that outside the SAA, as shown in Fig. [7J 

Each distribution is well explained by a Gaussian distribution. The fitting parameters of the 
Gaussian center and width are summarized in Table [IJ together with their errors. The ionizing 



Table 1: Result of Gaussian fitting to the histograms in Fig. [7J 



Center (ADU) Width (ADU) 

Outside SAA 151.70±0.04 13.26±0.04 
Inside SAA 175.90±0.05 13.82±0.05 
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radiation effects make the distribution shift toward higher intensities and increase its width (i.e. 
noise) significantly. But it keeps the symmetric shape with respect to the peak, which indicates 
that the rise of the averaged intensities inside the SAA is caused by a slight increase in intensities 
of all the pixels, but not an increase in the number of positive spiky pixels. We have confirmed the 
presence o f residual small posi t ive an d negative spikes below the detection limit of the de-spiking 



algorithm (jlshihara et al.ll2006l . I2010I ). producing the low-level broad wings in Fig. 7. However such 
spikes are only a few percent in number, which cannot make a significant contribution to the 
observed peak shift of the distribution. 

In Fig. [8l we plot the Gaussian widths against the central intensities, which are obtained by 
fitting a Gaussian profile to each histogram; by using the data taken between June 08 and June 14 
2006 (UT), we obtain many histgrams such as shown in Fig. If the offset of the output level is 
caused by a simple drift of the ROIC, the data taken inside the SAA should shift towards higher 
intensities along the horizontal axis while the Gaussian width does not change. However, the data 
taken inside the SAA show increases in the Gaussian widths as well as the intensities. Thus the 
noise does change due to the ionizing radiation effect in the SAA. Moreover the overall distributions 
of the data taken inside the SAA and outside the SAA show almost no systematic difference. Since 
dark current and photo-current are expected to have similar dependence of noise on intensities, 
Fig. indicates that the offset of an averaged pixel intensity inside the SAA is caused mainly by 
the actual increase in dark current, but not by a simple drift of the offset of the ROIC. 

Therefore, we conclude that the offset of the output level during the passage of the SAA is 
caused mainly by increases in dark current of the detector. 



5. Summary 

We have investigated ionizing radiation effects on the Si:As IBC mid-infrared detector on board 
the AKARI satellite. By using the AKARI all-sky-survey data, we have found that the detector is 
significantly affected by in-orbit ionizing radiation during the passage of the SAA. The effects are 
described as changes mostly in the offset of detector output, but not in the gain. We have confirmed 
that the distribution of the offset value shows an excellent spatial correlation with the distribution 
of the particle hitting rate in the terrestrial coordinate map, exhibiting no significant after-effects 
due to the passage of the SAA. The offset values show a linear correlation with the particle hitting 
rates. By utilizing the linear correlation, we have established a method to correct the ionizing 
radiation effects and successfully applied it to the all-sky-survey data, resulting in producing better 
diffuse emission maps. By investigating the histograms of pixel intensities in image units taken 
inside and outside the SAA, we conclude that the main physical cause of the offset variation is the 
increase of dark current. 



This research is based on observations with AKARI, a JAXA project with the participation 
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Fig. 2. — Difference of the foreground emission between two seasons, AI Z , plotted as a function 
of time during one revolution of the AKARI satellite for each of 5 sequential orbits. The triangles 
show observed AI Z , part of which are affected by the SAA or lost by pointed observations. The 
dashed lines show AI Z where the SAA-affected or lost data are estimated through interpolation by 
the observational AI Z for 24 hours. The perpendicular line roughly indicates the center of the SAA 
region. See text for details. 



- 11 - 




I i i i i i i i i I 

50 100 150 200 250 300 350 400 450 

Intensity (ADU) 

Fig. 3. — Fractional increase of the output level due to ionizing radiation effects {f{lf AA ) in eq.2) 
plotted against the averaged pixel intensities inside the SAA (lf AA ). We use the same datasets as 
those used in Fig. [2j The data points with the different ranges of particle hitting rates are indicated 
by different marks: the triangles for 40-60 unit -1 , circles for 80-100 unit -1 , squares for 120-140 
unit" 1 , and fitted by different hyperbolic curves. 
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Fig. 4. — Distribution map of the offset value of the detector output (I SAA — /) shown in the 
terrestrial coordinates. The contours show the distribution of the particle hitting rate, which are 
drawn by the 4 levels from 5 to 120 unit -1 on a linear scale. The lowest contour corresponds to 
the boundary of the SAA defined in this paper. The red line shows the trajectory of AKARI ihaX 
passes the SAA in the direction as indicated by the arrow. 
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Fig. 5. — Offset values of the averaged pixel intensities (I s AA — X) plotted as a function of the 
particle hitting rate. We use the same datasets as those used in Fig. [2j The dashed line shows a 
linear relation between the offset values and the particle hitting rates obtained by a line fittng. 
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Fig. 6. — AKARI 18 /im maps near the Galactic center: (a) original image created by including the 
data taken inside the SAA, (b) the image created after correcting the SAA-affected data by using 
the correction method presented in this paper. 
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Pixel intensity (ADU) 

Fig. 7. — Histograms of pixel intensities in the two image units viewing the same blank sky of 
10' x 48' area around (R.A., Dec.) = (159°. 8, — 25°. 0). The diamonds indicate the data taken 
inside the SAA, while the squares indicate the data taken outside the SAA. 
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Fig. 8. — Gaussian widths plotted against the Gaussian centers, both obtained by fitting to the 
histgrams of pixel intensities of image units. One data point is obtained from one histgram such 
as shown in Fig. [7J The data are taken between June 08 and June 14 2006 (UT). The diamonds 
indicate the data taken inside the SAA, while the squares indicate the data taken outside the SAA. 



